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Abstract 

The propagation of damage on the square Ising lattice with a corner geometry is studied by means of Monte 
Carlo simulations. By imposing free boundary conditions at which competing boundary magnetic fields ±/i 
are applied, the system undergoes a filling transition at a temperature Tf(h) lower than the Onsager critical 
temperature Tc. The competing fields cause the formation of two magnetic domains with opposite orientation 
of the magnetization, separated by an interface that for T larger than Tf(h) (but T < To) runs along the 
diagonal of the sample that connects the corners where the magnetic fields of different orientation meet. Also, 
for T < Tf(h) this interface is localized either close to the corner where the magnetic field is positive or close to 
the opposite one, with the same probability. 

It is found that, just at T = Tf(h), the damage initially propagates along the interface of the competing 
domains, according to a power law given by D(t) oc t v . 

The value obtained for the dynamic exponent (rf = 0.89(1)) is in agreement with that corresponding to the 
wetting transition in the slit geometry (Abraham Model) given by 77 T = 0.91(1). However, for later times 
the propagation crosses to a new regime such as rj** = 0.40 ± 0.02, which is due to the propagation of the 
damage into the bulk of the magnetic domains. This result can be understood due to the constraints imposed 
to the propagation of damage by the corner geometry of the system that cause healing at the corners where the 
interface is attached. 

The critical points for the damage spreading transition (To (ft)) are evaluated by extrapolation to the ther- 
modynamic limit by using a finite-size scaling approach. Considering error bars, an overlap between the filling 
and the damage spreading transitions is found, such that Tf{h) = To(h). 

The probability distribution of the damage average position P(lo) and that of the interface between magnetic 
domains of different orientation P(lo) are evaluated and compared. It is found that, within the nonwet phase, 
the average position of the damage lies slightly shifted from the interface toward the side of the largest domain. 
However, in the wet phase both P(l ) and P(lo) are Gaussians exhibiting a single peak at the position of the 
diagonal of the corner sample. 
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1 Introduction 



The interplay between critical behavior and confinement is a challenging topic in the field of Condensed Matter 
and Statistical Physics. In fact, the confinement of fluids, polymers, magnetic materials, etc., by walls that interact 
with the physical system, leads to the occurrence of very interesting phenomena such as capillary condensation, 
wetting, corner wetting (filling), etc. PQG]- 

Particularly interesting is the case of confinement in two dimensions (d = 2) close to second-order critical 
phenomena, due to the occurrence of strong fluctuations. Within this context the classical Ising Model has been 
taken as an archetypical system for the study of capillary condensation, wetting and corner filling by means of 
analytical techniques (capillary wave theories, transfer matrix method, etc.) and extensive numerical simulations 
[HI HI HI EH IHj (for a review see e.g. [2| and references therein). 

Wetting is usually studied by using the strip geometry shown in Figure^] In this case, the presence of competing 
magnetic fields (h) along the confinement walls induces the formation of an interface between domains where most 
spins are pointing up and down. This interface runs essentially parallel to the walls, and in a finite system, it 
undergoes a localization-delocalization "transition" that is the precursor of a true second-order wetting transition 
that takes place in the thermodynamic limit only (L — ► oo, M — > oo). The phase diagram (i.e. the critical curve in 
the h — T plane, as shown in Figure has been solved exactly by Abraham [HIE], yielding 

cosh(2/i/3) = cosh(2if ) - e- 2K sinh(2iT), (1) 

where J > is the coupling constant, h is the surface magnetic field, j3 = 1/kT is the Boltzmann factor, and 
K = J/3. 

On the other hand, the study of corner filling is performed by using the geometry sketched in Figure El where 
the competing magnetic fields are applied at opposite corners. The study of this filling transition under equilibrium 
conditions has recently attracted growing attention fTTfl fill ril^rHfl fill rHflrM 

ESI EOl EU E21 E31 EH EH ESI E3 EB] Also, the filling transition upon the irreversible growth of a magnetic system 
has very recently been studied In both cases, the occurrence of an interface is due to the presence of competing 
fields. The localization-delocalization transition of the interface in a finite system yields to a true second-order 
corner filling transition in the thermodynamic limit (L — > oo). The analytical expression of the equilibrium phase 
diagram was early conjectured by Parry et al [22] and more recently proved rigorously by Abraham and Maciolek 
[23, yielding 

cosh(2/i/3) = cosh(2iT) - e~ 2K sinh 2 (2K). (2) 

The critical curve is presented in Figure showing that for a given surface magnetic field, the filling transition 
takes place at a lower temperature than the wetting transition, except of course for h = 0, where both curves 
converge to the Onsager critical temperature of the Ising Model (Tc). 

In previous work we showed that the presence of interfaces between magnetic domains of different orientations, 
such as those observed close to the wetting transition, favors the propagation of perturbations in magnetic materials 
[30IE1E21- A standard technique frequently used to study the propagation of perturbations is the measurement 
of the Hamming distance or Damage (D(t)), given by 

1 N 

£(*) = ^£|^(i,r)-Sf(t,T)|, ( 3 ) 
/ 

where the summation runs over the total number of spins TV, and the index I (1 < I < N) is the label that identifies 
the spins of the configurations. S (t, T) is an equilibrium configuration of the system at temperature T and time 
f, while S B (t,T) is the perturbed configuration that is obtained from the previous one just by flipping few spins 
[331 B3j. In order to further contribute to the understanding of the propagation of perturbations in magnetic 
materials and to clarify the role played by both the interfaces and the confinement geometry, the aim of this work 
is to report an extensive computer simulation study of damage spreading close to the corner filling transition of the 
d = 2 Ising Model. Furthermore, the obtained results will be discussed within the context of our current knowledge 
of the corner filling transition [23 and compared with our previous studies of damage propagation in the strip 
geometry where a wetting transition is observed |401 14H 142] . While in both cases one has fluctuating interfaces 
due to the presence of competing fields, the geometries used are quite different (see e.g. Figures and for the 
sake of comparison) , so we expect that this situation will help us to contribute to the understanding of the effect of 
confinement on the properties of interfaces and the propagation of damage along them. 
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This manuscript is organized as follows: in Section II we describe the corner wetting transition and the Monte 
Carlo simulation method. In Section III we present and discuss the obtained results. Finally, the conclusions are 
stated in Section IV. 



2 Brief overview of the corner filling transition and the simulation method 

The Hamiltonian (TL) for the Ising Model in the corner geometry with competing short-range local fields at the 
boundaries, as sketched in Figure is given by 

TL = -J Vi,j(7m,n - - hy^OLj + h } j (Tij + h}a itL , (4) 

<i, i 7,m,n> i j j i 

where aij = ±1 is the spin variable, J > is the coupling constant, and h is the magnitude of the surface field. 
The first summation runs over all spins, while the remaining ones hold for spins at the surfaces where the magnetic 
fields are applied (see also Figure [2j and h > is measured in units of J. 

It is very useful to discuss the filling transition of a cavity by using the mapping between the Ising Model and 
a Lattice Gas. So, by considering the transformation spin up ^ vapor, spin down liquid, respectively, one has 
that the derealization of the interface can be thought as the growth of a macroscopic liquid layer wetting on the 
substratum. Let us consider a corner cavity forming an angle ±<fi with the horizontal axis, as shown in Figure 
|2j We further assume that the system is in contact with the vapor phase of the bulk at a certain temperature T 
and chemical potential fx, so that at coexistence one has fi = (i S at- According to (macroscopic) thermodynamic 
arguments it has been shown ^J] that complete filling by the liquid is observed for T >Tf, where Tf is the filling 
temperature given by 

9(T f ) = 0, (5) 

where 9(Tf) is the contact angle of the liquid to the flat substratum. Consequently, one has that Tf < T w , where 
T w is the critical temperature for the wetting of the plane by the liquid. 

The case studied here corresponds to <f) — n/4 (see Figure EJ- The corner filling transition of the Ising Model 
is of second order, so one has power-law behavior of relevant quantities such that the average interface position in 
the y-direction ((Iq)), the parallel and perpendicular (£j_) correlation lengths (see Figure that describe the 
fluctuations of the interface in their respective directions, behave as follows 

(Jo) oc e- ft ; £l cx e -"± ; oc e""-, (6) 

where e = {hc{T) — h) is the distance to the critical curve, /3 S = 1 is the order parameter critical exponent |45j . 
and v± = v x = 1 are the correlation length exponents in the direction perpendicular and parallel to the interface, 
respectively |29] , 

On the other hand, the probability distribution of the position of the interface is given by 

P^oc^f^-^Y (7) 



(lo) \(lo)'{lo) 

where the second scaling variable can be neglected close to criticality because both (l ) and £j_ diverge with the 
same critical exponent (see equation ©). Furthermore, at criticality and for L — > 00, e — > one has ^2] 

P(fo)oc-^e X pHo/(Zo», (8) 
(to) 

indicating that the probability of finding the interface close to the walls decays exponentially. For large enough 
lattices, equation © can be symmetrized so that 



PQo) oc ' 



2(/o> 



exp —— + exp 1 



.(h) J ' ' V (/o> 

All these equations hold for the localized interface below but close to the critical curve shown in Figure 
On the other hand, at the critical point the distribution is expected to be flat ^2], namely 



(9) 
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P(Z ;e = 0) = i, L^oo, (10) 

so that the interface can be located at any place with the same probability. 

Above the filling transition the average position of the interface lies along the diagonal of the sample with 
(la) = L/2 [46.', and then the distribution is a Gaussian given by 

P(*o)ocexp(-fc^) (11) 

with £j_ oc . For further details on the corner filling transition see and references therein. 

The spreading of damage in the corner geometry is studied by means of Monte Carlo simulations assuming the 
Glauber dynamics. So, a randomly selected spin is flipped with probability p(flip) given by 02] 

exp(-/3 ■ Aff) 

Pifhp) = l + exp(-/3-Ag) ' (12) 

where AH is the change of the Hamiltonian given by equation <@J due to the attempted flip and (3=1/ kT is the 
Boltzmann factor. In order to set the time scale, we assume that during a Monte Carlo time step (mcs) all spins of 
the system (L x L in total) have the chance to be flipped once, on the average. 



3 Results and Discussion 

Simulations are performed on the square lattice of size L x L (64 < L < 2048) by using the geometry shown in 
Figure Eland by assuming open boundary conditions. Most simulations are performed by keeping the temperature 
constant (0.75 < T/Tc < 0.95) and changing the magnitude of the surface field (0.1 < h < 0.35). These ranges 
for the parameters are selected to avoid corrections due to the standard critical behavior of the Ising magnet that 
appear for T — > Tq and h — ► 0. 

The dynamic behavior of the system is characterized by measuring the total damage or Hamming distance 
according to equation <j2j and the probability distribution of the distance between the damage interface and the 
corner P(Iq) (see Figure EJ. 

In order to obtain equilibrated configurations, the ground-state configuration corresponding to T = (all spins 
pointing up) was annealed up to the desired set (h,T) for 10 4 mcs. Subsequently, initial damage was created 
by flipping the spins of the diagonal perpendicular to the interface (i.e. along the vertical direction in Figure El 
according to this general rule: if the magnetization of the whole sample is positive (negative), the up (down) spins 
are flipped. Using this procedure, it is assured that the initial damage is always D(t = 0) < L/2 ^fjj. This kind 
of perturbation reproduces the effect of a large magnetic field applied at the diagonal of the sample and pointing 
away to the direction opposite to the total magnetization. 

Subsequently, the time evolution of the damage was recorded. FigureQ]shows plots of D(t) versus t obtained by 
keeping T = 0.80Tc constant and varying the surface field h. The observed behavior resembles the results, already 
published, corresponding to the wetting transition in a strip geometry j321- In fact, the short-time behavior of the 
curves involves the healing of the damage initially created in the bulk of the domains that takes place up to t ~ 50 
mcs. Subsequently, the propagation of the damage along the interface becomes relevant. While for weak fields one 
observes damage healing (e.g. h < 0.35), by increasing the field the damage starts to propagate monotonically and 
a power-law behavior of the form 

D(t)<xt v , (13) 

where T) is an exponent, can be proposed as in the case of the strip geometry [32] . 

The onset of a power-law divergence, for a sample of size L, is considered as the signature for a size-dependent 
"critical" behavior of the damage process. In fact, as shown in Figure the power-law behavior is observed at 
slightly different surface magnetic fields when the size of the lattice is changed. Moreover, simulations performed 
using lattices of different size confirm that the power-law behavior of the damage becomes more evident for larger 
samples, as shown in Figure HJ These results anticipate that one would need to perform a proper extrapolation of 
the data in order to obtain the actual critical points in the thermodynamic limit. Before tackling this issue, let us 
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first discuss the behavior of D(t) as shown in Figure EJ An overview inspection reveals that three different time 
regimes can easily be distinguished, as follows: 

1) During the short-time regime (t < t m i n ), the damage decreases monotonically reaching an absolute minimum 
at t m in, which only depends slightly on the lattice size (see dashed-dotted line in Figure El- This behavior can be 
understood by recalling that the damage is initially generated along a diagonal line of the sample, in the direction 
perpendicular to the spin up-spin down interface (see Figure GJ|. This initial condition creates - almost linear - 
damaged regions in the bulk of well-ordered magnetic domains (see the snapshot of Figure taken for t — lmcs) . 
Due to the very low density of broken bonds in the bulk, the damage cannot propagate and becomes quickly healed 
(see the snapshot of Figure taken for t — 5mcs) while a small fraction of the initially created damage still remains 
at the interface (see the snapshot of Figured taken for t = 28mcs). Since the initial damage is distributed along the 
line one has D(t = 0) ~ in agreement with the equispaced intercepts shown in FigureEJfor t = 0. However, at 
t m in one has only few damaged sites along the center of the interface, so one expects D(t m in) ~ L ~ 2 , a result that 
is corroborated by the numerical data as shown in the inset of Figure Summing up, the initial decrease of the 
damage is due to damage healing in the bulk. 

After reaching t m in one observes the growth of the damage. This behavior can be better observed in the scaled 
plot of Figure that attempts to collapse all curves corresponding to different lattices. Recall that a perfect 
collapse cannot be achieved not only due to the existence of prefactors but also because D(t = 0) ~ L^ 1 while 
D(t m in) ~ L~ 2 . However, the scaling plot of Figure^lis useful in order to clearly show the power-law divergence of 
the damage above t m in- Here, two well-defined regions can also be distinguished (recall that these regions become 
also evident in Figure EJ. 

2) For t > tmin the damage increases according to equation ((13)) and the best fit for the data gives rf — 0.89(1) 
(see the dashed line in Figure El- It is worth mentioning that this exponent is in good agreement with our previous 
measurement of damage spreading along the interface in the strip geometry (see Figure^ at the wetting transition, 
which yields rj WT — 0.91(1) 021- So, we conclude that the critical increase of the damage, which takes place for 
tmin < t < t cross , where t cross is the crossover time to a third regime that will be discussed below, is due to the 
propagation of damage along the spin up-spin down interface characteristic of the filling transition. This statement is 
further confirmed by the snapshot configuration shown in Figure0for t = 500mcs. Then, one has that the exponent 
describing the propagation of the damage along the interface between magnetic domains of opposite magnetization 
is given by r/ 1 = 0.90(2), where the error bars are large enough to account for the fact that rj 1 — r/ WT = -q* . Our 
findings are also consistent with the fact that the propagation of damage along the interface may be independent 
of the nature of the studied phenomena - wetting or filling - provided that the measurements are performed at 
criticality. 

3) Let us now discuss the propagation of damage after the crossover time t cross . Here, a power law is also observed 
(see the dotted line in Figure ^Jl but the best fit of the data according to equation 11 -til yields rj** — 0.40(2). This 
result, measured at the critical point, depends neither on L nor on hn(L), So, we conclude that this exponent 
characterizes the propagation of the damage from the interface into the bulk of the corner geometry (see the snapshot 
of Figure taken for t = 50000mcs) . At this latter stage the propagation of the damage along the interface has 
ceased due to the constraint imposed by the corner geometry and one observes a slower propagation (rj 1 > r/**) into 
the bulk adjacent to the interface. Also notice that the damage generated along the interface and close to it follows 
the fluctuations of the actual position of the interface. 

Of course, during the second time regime (t m in < t < t cross ) the damage does not strictly propagate along the 
interface only, since one also expects the onset of propagation in the direction perpendicular to it. However, due to 
the fact that the magnitudes of the exponents are quite different, this effect is no longer observed in the simulations. 

After discussing the dynamics of the damage, we will like to focus our attention on the location of the damage 
spreading critical points in order to draw the corresponding phase diagram. For this purpose, note that the best 
fit of the power-law behavior according to equation il'M) - within the long-time behavior for t > t cross - allows us 
to identify the size-dependent critical field (/jd(L)) for a given temperature (see e.g. FigureEJthat corresponds to 
T = 0.80Tc). It is well known that in numerical simulations, performed by using finite samples, "critical" points 
are shifted and rounded due to finite size effects. So, as in the case of the wetting transition jH], we propose the 
following Ansatz for the shifting of the damage-healing critical point 

ho(L) - h D (oo) ~ L~ 7 , T=constant (14) 
where 7 is an exponent and ho (00) is the true critical field in the thermodynamic limit. 
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The best fits for the data corresponding to four different temperatures, as shown in Figure 03 were obtained by 
taking 7 = 2. The extrapolated values of the critical field are listed in Table 1 that also includes, for the sake of 
comparison, the critical fields obtained by solving the exact solution of the corner-filling phase diagram 00] given 
by equation 



T/T c 


h D (oo) eq. Ili4l) 


hf(oo) eq. (UJ| 


0.75 


0.400 ±0.005 


0.3975 


0.80 


0.360 ±0.005 


0.3505 


0.90 


0.245 ±0.005 


0.2411 


0.95 


0.175 ±0.005 


0.1683 



Table I. Critical field for damage spreading (/i_d(oo)) obtained by extrapolation to the thermodynamic limit 
with the aid of equation 1)1411 . The third column corresponds to exact values of the critical field correspondig to the 
filling transition as obtained by using equation (|2J). 

As follows from Table 1 and Figure the critical points for damage spreading in the corner geometry are 
indistinguishable from those of the exact solution of the filling transition (within error bars). This result is in 
contrast with our previous studies close to the wetting transition in a strip geometry where one has hn(oo) < h w (oo) 
(see also Figure and the damage spreading transition is located in the nonwet phase of the phase diagram. We 
expect that this behavior may be due to geometrical constraints imposed by the corner array. In fact, in the strip 
geometry and for h < h w (co) one has that the interface is still bound to one of the walls. However, even in this 
nonwet phase but close to the wetting transition, the damage propagates along the interface without geometrical 
constraints. On the other hand, in the corner geometry when the interface is bound to one corner, the spatial 
propagation of the damage is restricted. 

In order to gain further insight into the spatiotemporal propagation of the damage we also measured the 
probability distribution of the distance from the damage zone to the corner (P(Iq)). The distribution was evaluated 
along the diagonal of the sample as shown in Figure Figure El shows a summary of the obtained results. 

For h <C hf (00) (e.g. h = 0.20) the distribution is almost flat with two small peaks close to the corners. 
Approaching the transition by increasing the field these peaks develop and become slightly shifted toward the 
center of the sample (e.g. h = 0.21 and h — 0.22 in Figure EJ). This double-peaked structure indicates that the 
damage remains bound to each corner with the same probability as expected for the case of the nonwet phase. On 
the other hand, for h ~ hf (00) (e.g. h = 0.23, 0.24 in Figure^ the distribution becomes a Gaussian centered along 
the middle of the sample. The Gaussian structure of P(Iq) remains even for h 3> hf (00) (e.g. h = 0.30 in Figure 

It is also very useful to compare P{Iq) with the probability distribution of the position of the interface between 
domains of opposite magnetization P(lo). Figure llOl shows plots of both P(Iq) and P(^o) obtained far below the 
filling transition. One observes that for the domain interface P(lo) decays exponentially, as expected according 
to equation ©, indicating that for this set (h,T) such an interface is strongly localized close to each corner with 
the same probability. On the other hand, while the interface of the damage is still bound, its average position lies 
further apart from each corner, as evidenced by the double-peaked structure of P{1$) (see Figure El- So, Figure 
ITol provides clear evidence that the damage is located in the neighborhood of the interface between domains but 
slightly shifted toward the bulk of the sample, or more specifically within the largest domain. 

Close to the filling transition (FigureEJ ° ne h as that the interface between domains exhibits a flat distribution, 
as expected from equation (1101 . indicating that the interface could be found at any place of the sample. However, 
P(Iq) exhibits a single peak centered at the middle of the sample reflecting the inertia of the damage in order to 
follow the displacement of the interface. Furthermore, the abrupt decay of the damage close to the corners, which 
is evidenced in both Figures ITfil and ITT1 is due to the healing effect caused by the neighboring magnetic fields. This 
effect prevents the damage from approaching the corners acting as an effective repulsive effect. Also, this result is 
consistent with our previous conclusion that the damage is essentially located within the largest magnetic domain. 

Finally, within the wet phase both P(lo) and P(lo) are centered around the diagonal of the sample (see Figure 
IT2I1 in agreement with the fact that the interface is delocalized . 

The main conclusions that follow from the comparison between P{Iq) and P{lo) (see FiguresElElEl are, on 
the one hand, the operation of an effective repulsion of the damage at the corners where the fields have the same sign, 
and, on the other hand, that the damage becomes attached to the interface but preferentially located toward the 
largest domain. These results are in contrast to previous observations performed by studying the spreading of the 
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damage close to the wetting transition using the strip geometry where the damage was uniformly and symmetrically 
distributed along the interface between competing domains 02] • We expect that this difference between corner and 
strip geometries may explain the fact that for the former the phase diagram for damage spreading matches that of 
corner filling (Figure while for the latter the damage critical points lie within the nonwet phase (Figure El- 
Let us now recall that, within the wet phase, the Gaussian distribution of the interface profile (equation ltTT|l l 
translates into a Gaussian distribution of the magnetization m given by [22j 

P L (m) oc exp( ) (15) 

where \l is the magnetic susceptibility and (3 — 1/kT is the Boltzmann factor. Equation 11511 simply reflects the 
fact that the interface is located, on average, along the diagonal of the sample, the magnetic domains of different 
magnetization having the same average size. On the other hand, the susceptibility diverges with the lattice size 
according to xl °c (M 2 ) oc L In view of these facts we also measured the width of the Gaussian distributions 
P(lg) (see the full line in Fi gure I12|l within the wet phase, given by {sfj}. The inset of Figure El shows log-log 
plots of (sfj) versus [h — hf] obtained by keeping T — 0.907c constant and for different lattice sizes. A preliminary 
inspection shows that (sf-,) increases when the field approaches the critical value, as well as when larger lattices 
are considered, suggesting that (sp) ~ xl, within the wet phase. In order to test this observation, it is worth 
mentioning that starting from the distribution of the magnetization (equation @H ) it is possible to obtain the 
standard relationship between the susceptibility and the fluctuations of the order parameter given by (22] 

k B T X = L 2 ((m 2 ) - (H) 2 ) - L 2 ■ x{Le), (16) 

where \m\ is the absolute value of the magnetization, x{z) is a scaling function, and e = (hf — h). On the other 
hand, within the wet phase, one has e < and (\m\) = (at least in the thermodynamic limit). Then, if the 
proposed proportionality ((s|>) ~ xl) holds, one would have 

(s 2 D ) ~ L 2 ■ f(Ls), (17) 

where / is a scaling function. Figure shows the corresponding scaling plot (s 2 D )L~ 2 versus L{h — hf)) as obtained 
by using the data of the main panel. By considering the involved errors and the relatively small samples that can be 
simulated with our computational resources, we conclude that the collapse is acceptable and most likely equation 
1(17)1 should hold. 



4 Conclusions 

We studied the spreading of damage in the two-dimensional Ising Model confined in a corner geometry with a 
competing (short-range) magnetic field acting on the surfaces (see FigureGJi- By considering initial damage, created 
along one diagonal of the sample, in the direction perpendicular to the interface between domains of opposite 
magnetization originated by the applied fields, we conclude that the dynamics of damage spreading exhibits three 
characteristic regimes: i) For short times, one observes the healing of the damage created in the bulk of the 
domains. However, at the critical damage spreading point, a small cluster of damaged sites always survives close to 
the interface, ii) This small cluster propagates along the interface according to a power-law behavior D(t) oc t n , 
with rj* = 0.89(1). By comparing the results of the present work and those already published 02] for the spreading 
of damage along the interface at the wetting transition, we conclude that the exponent r] 1 = 0.90(2) describes the 
critical propagation of the damage along magnetic interfaces, iii) Finally, due to the constraint imposed by the 
corners where magnetic fields of opposite direction meet, the damage no longer propagates along the surface but 
starts to slowly spread into the bulk of the domains. Within this regime one has D(t) oc f *, where 77** = 0.40(2) 
is the exponent describing the spreading of the damage in the bulk. 

We located the damage spreading transition by proper extrapolation to the thermodynamic limit. It is found 
that, within error bars, the phase diagram coincides with that of the corner filling transition. This result is in 
contrast to the previous study of the wetting transition where the damage spreading transition lies within the wet 
phase. This difference can be understood on the basis of the constraint imposed by the field acting at the corners, 
causing a repulsive effect on the damage. This repulsive effect is nicely observed in measurements of the probability 
distribution of the damaged area within the nonwet phase. Finally, in the wet phase the distribution of the damage 
is Gaussian and ist width scales as the fluctuations of the magnetization. 
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We hope that this study will contribute to the understanding of the propagation of magnetic perturbations along 
magnetic interfaces, as well as to the understanding of more complex damage spreading transitions that depart from 
the universality class of directed percolation. 
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Figure 1: Strip geometry of size L x M. The surface magnetic fields of different orientations applied at the upper 
(bottom) rows of the lattice are indicated by + (— ) signs. Notice that periodic (open) boundary conditions are 
assumed along the x (y) axis. See notation and further details in the text. 
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Figure 2: Comer geometry of size L x L. The signs + and — indicate the surfaces where the competing surface 
magnetic fields are applied. In this case, the boundary conditions are open for all sides of the sample. See notation 
and further details in the text. 
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Figure 3: Phase diagram of h versus T/Tc- The dashed line corresponds to the wetting transition (equation dJ) 
while the full line corresponds to the corner filling transition (equation @). The circles are obtained for damage 
transition at the strip geometry (32 an d are shown for the sake of comparison. The squares correspond to the 
results obtained in the present work for the damage spreading transition at the corner geometry. 
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Figure 4: Log-log plot of damage (D(t)) as a function of time. Results obtained for different values of the surface 
magnetic fields h, as listed in the figure. The simulations were performed by using a lattice of size L = 256 and for 
T = 0.80T C . The dashed line has slope 77** = 0.40. 
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Figure 5: Log-log plot of D(t) versus t, obtained for T = 0.80Tc and at the size-dependent "critical" magnetic field 
h D (L): h D (L = 256) = 0.355, h D (L = 512) = 0.3675, h D (L = 1024) = 0.37, and h D (L = 2048) = 0.38. The 
dashed line has slope rj* = 0.89 and the dotted line has slope rj** — 0.40. The dashed-dotted line shows the slight 
dependence of the minimun value of the damage (D(t m i n )) on L. The inset shows a linear-linear plot of D(t m i n ) 
versus L~ 2 . 
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Figure 6: Scaling plot (D(t) ■ L 2 versus t) of the data shown in Figure El The departure of the data observed for 
t — ^ is consistent with the initial condition involving the generation of a damaged line with D(0) ~ L~ . The 
dashed line has slope rj* = 0.89 and the dotted line has slope 77** = 0.40. 
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Figure 7: (a) Snapshots of equilibrium configurations obtained at T — 0.8Tc, h = 0.37, L = 256, and different times 
from top to bottom: t = 1 mcs, t = 5 mcs, t — 28 mcs, i = 500 mcs, and t = 50000 mcs. (b) Snapshot pictures of 
the damaged sites obtained from the lattices shown in (a) by applying linear damage along the diagonal. Damaged 
sites are shown in black and undamaged sites are left white. 
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Figure 8: Plots of the size-dependent "critical" magnetic fields of the damage transition hn{L) as a function of L 2 , 
as obtained for the temperatures listed in the Figure. 
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Figure 9: Plot of probability distribution of the position of the damage P(Iq), obtained for L = 128, T = 0.901c, 
and different values of surface magnetic field h, as listed in the figure. 
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Figure 10: Probability distribution of the interface position P(lo) (full line) and damage P(Iq) (dotted line) along 
the y-direction. Data corresponding to L — 128, T = 0.90Tc, and h = 0.22 < hf(oo). 
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Figure 11: Probability distribution of the interface position P(lo) (full line) and damage P(lo) (dotted line) along 
the y-direction, for L = 64, T = 0.90T C , and h = 0.24 w h f (oo). 
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Figure 12: Probability distribution of the interface position P(7o) (squares) and damage P(Iq) (circles) along the 
y-direction, for L — 128, T = 0.90Tc, and h = 0.25 > /i/(oo). The full line corresponds to the Gaussian fit of the 
data. The data corresponding to P(lo) were multiplied by a factor of 5 for the sake of clarity. 
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Figure 13: Width of the distribution of damage (P(Iq)) (s%) versus [h — hf], obtained at T = 0.90Tc and within 
the wet phase. Results corresponding to different values of the lattice size L, as indicated in the Figure. Inset: 
Collapse of the data according to equation II17II . 
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